Re-examination of Polytropic Spheres in Palatini f(R) Gravity 
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We investigate spherically symmetric, static matter configurations with polytropic equation of 
state for a class of f(R) models in Palatini formalism and show that the surface singularities recently 
reported in the literature are not physical in the case of Planck scale modified lagrangians. In such 
cases, they are just an artifact of the idealized equation of state used. In fact, we show that for 
the models f(R) = R ± XR 2 , with A on the order of the Planck length squared, the presence of a 
single electron in the Universe would be enough to cure all stellar singularities of this type. From 
our analysis it also follows that the stellar structure derived from these lagrangians is virtually 
undistinguishable from that corresponding to General Relativity. For ultraviolet corrected models 
far from the Planck scale, however, the surface singularities may indeed arise in the region of validity 
of the polytropic equation of state. This fact can be used to place constraints on the parameters of 
particular models. 

PACS numbers: 98.80.Es , 04.50.+h, 04.25.Nx 
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I. INTRODUCTION 



The study of modified theories of gravity is interesting 
for many reasons. General Relativity (GR) has stood as 
a monument to the genius of Albert Einstein for almost a 
century. However, despite the astonishingobservational 
and experimental success of this theory there exist 
many theoretical aspects which make us consider the 
possibility of departures from it in certain regimes. In 
fact, one would like to remove physically undesired 
aspects such as the Big Bang or black hole singularities. 
These are aspects which clearly point towards high en- 
ergy modifications of the theory motivated by quantum 
gravity effects. On the other hand, the recently observed 
cosmic speedup @ suggests that something is needed in 
the Friedman-Robertson- Walker equations to account 
for the apparent positive acceleration of the expansion 
rate. This requires a modification of the equations 
at low energies. It could come in the form of a new 
energy source dominant over the other matter sources 
at very low densities @ , or as a modification of the field 
equations which would be suppressed above a certain 
cosmic curvature scale Q. It is this latter possibility 
which has impulsed the study of modified theories of 
gravity of the f(R) type in the recent literature [f| (see 
also (1). 

When the Einstein-Hilbert lagrangian is promoted to a 
non-linear function of the scalar curvature one is dealing 
with an f(R) theory of gravity. When this is done, 
one must specify whether the connection is still defined 
in terms of the metric as the Levi-Civita connection 
(metric formalism) or if it is an independent geometrical 
field (Palatini formalism). One should note that, 



on pure geometrical grounds, metric and connection 
are completely independent and equally fundamental 
objects. However, the literature has traditionally given 
a dominant role to the metric over the connection. In 
this paper we will focus on f(R) theories in Palatini 
formalism. In these theories, the field equations show 
that the independent connection is not a dynamical 
field in the sense that it does not satisfy an independent 
second-order partial differential equation. Instead, its 
field equation can be solved in terms of an auxiliary 
metric which is conformally related with the usual (or 
physical or Jordan frame) metric. As we will see, the 
conformal factor is a local function of the trace, T, of the 
energy-momentum tensor T^ v of the matter. Thus, the 
independent connection, despite being non-dynamical in 
the usual sense, does mediate new interactions between 
the matter and geometry. The resulting equations 
for the metric can be written in Einstein form with a 
modified right hand side in which the trace T plays a role 
non-existing in GR. We will discuss later in depth the im- 
plications of such terms in the dynamics of these theories. 

The goal of this paper is to contribute and shed some 
new light on a problem discussed in recent literature 
that seems to affect seriously the theoretical viability 
of all f(R) models in Palatini formalism. The problem 
consists on the existence of curvature singularities near 
the surface of spherically symmetric, static matter 
configurations with polytropic equation of state (EoS) 
and polytropic index within the range 3/2 < 7 < 2. Its 
existence and possible consequences were first reported 
by Barausse et al. in Q using previous results of 
Kainulainen et al. ||. The problem was reconsidered 
by Kainulainen et al. in Q, where it was argued that 
the singularities had more to do with the peculiarities 
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of the EoS used than with the theory of Palatini f(R) 
gravity. They found that for neutron stars the tidal 
acceleration due to the singularity becomes equal to the 
Schwarzschild value of GR only at a distance ~ 0.3 fermi 
from the surface of the star, which makes unrealistic the 
use of a polytropic EoS. However, Barausse et al. replied 
in [Toj that the fluid approximation is still valid on the 
scale at which the tidal forces diverge just below the 
surface of a polytropic sphere in the case of the generic 
functions f(R) considered. 

We reconsider the problem and independently re-derive 
the structure equations for spherically symmetric, static 
systems. We study the physically admissible models 
f(R) = R ± XR 2 , with A a constant of order the Planck 
length squared, and estimate the density scale p s at 
which the potentially divergent terms become of order 
unity. We find that, for our choice of A, p a is well below 
any physically accessible density scale. In this case, one 
can conclude that the mathematical singularity found 
in Q and is just an artifact due to the use of the 
polytropic EoS beyond its natural regime of validity. 
For other choices of A, however, the singularity can 
effectively arise at densities sufficiently high as to admit 
the validity of the polytropic EoS. Our approach and 
results can thus be used to place bounds on the allowed 
range of values of A and to constrain other f(R) models. 

The paper is organized as follows. In section II we 
present the action, derive the field equations, and discuss 
their physical properties and new aspects as compared 
to GR. Section III presents the structure equations for 
a spherically symmetric, static matter configuration for 
a generic, unspecified lagrangian f(R). In section IV we 
study the behavior of the metric near the surface for sys- 
tems with polytropic EoS. We conclude with a summary 
and discussion. 



II. ACTION AND FIELD EQUATIONS 

Let us begin by defining the action of Palatini theories 

%, T, tp m ] = J di xy/^gf(R) + S m [g^ v , i/) m ] (1) 

Here f(R) is a function of R = g^R^T), with R^(T) 
given by R^(T) = -0„I*„ + d x T^ + T^T^ - r*,r£ A 
where T^ v is the connection. The matter action S m de- 
pends on the matter fields ip m , the metric g^, which 
defines the line element ds 2 = g^dx^dx" ', and its first 
derivatives (Christoffcl symbols). The matter action does 
not depend on the connection , which is seen as an 
independent field appearing only in the gravitational ac- 
tion (this condition is not essential and can be relaxed at 
the cost of introducing a non- vanishing torsion). Varying 
([!]) with respect to the metric g^ we obtain 



where f'(R) = df fdR. From this equation we see that 
the scalar R can be solved as an algebraic function of the 
trace T. This follows from the trace of |J2) 



f(R)R-2f(R) = K 2 T, 



(3) 



The solution to this algebraic equation will be denoted 
by R = TZ(T). The variation of {TJ with respect to T* u 
must vanish independently of @ and gives 



(4) 



where is the derivative operator associated to T and 
/' = f'(lZ[T]) is a function of the matter fields. This 
equation leads to 

a t Xp 

^fiu = {d^tpv + d u t pfl — dpt^) (5) 

where t^ v = fg^v, and /' = Jjj- Inserting this solution 
for r into Q we find the following modified Einstein 
equations 



If- 



(0) 



where now is the usual covariant derivative of gp V . 



A. Physical aspects of the field equations 

First of all, it is important to note that the right hand 
side of ® behaves like a modified energy-momentum 
tensor in which the trace T plays a role non-existing in 
GR via the terms 1Z, f(TZ) and the various derivatives of 
f(TZ). Thus, though the equations for the metric are still 
second order, the dynamics in regions with pronounced 
gradients in the matter fields (df and V M V„/') might 
considerably depart from that of standard GR. In 
contrast, the vacuum solutions of these equations are 
exactly the same as in vacuum GR with a cosmological 
constant. 

To better understand the role of the gradients of /' in 
©, it is useful to re-express the field equations in terms 
of the auxiliary metric that defines the connection T 
in ©. In terms of this metric, the equations become 



G>i/(t) — T^ v 



W - f] 

2(f) 2 ' 



(7) 



f(R)R^(T) - -f(R)gp V = n 2 T, 



(2) 



The advantage of using is that the gradients d^f 
and V^W/' present in ((6]) are completely absent in 
([7]), which allows for a better comparison with GR. 
The structure of these equations resembles that of GR 
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plus a trace-dependent term 



2(/') 



which becomes a 



cosmological constant in vacuum. The solution for t^ u 
will thus involve integrals over the matter sources which, 
like in GR, represent the fact that larger amounts of 
matter contribute more to the curvature of space-time. 
What is then the role of the gradients d^f and V^Vj,/' 
of ([6])? The answer is simple: since the physical metric 
is g^v = -ff > the gradients of /' in © introduce a local 
dependence of the metric on the local matter-energy 
distribution. The metric is not only determined by the 
total amount of matter, also the local matter-energy 
density affects its properties via f'(lZ(T)). This depen- 
dence on the local energy density distribution manifests 
itself with clarity in f(R) models sensitive to low energy 
scales. In particular, it has been shown that all such 
models would have dramatic effects for the stability 
of matter (Tl| . For the model f(R) = R - p, 4 /R, 
the ground state of Hydrogen would desintegrate in 
less than two hours [12|], which rules out that model 
and all models sensitive to ultra-low matter density 
scales. The reason for this has to do with the strong 
gravitational backreaction induced by /' near the zeros 
of the atomic wave functions, where the characteristic 
ultra-low density scale p M = p 2 / n 2 ~ 10~ 26 g/cm 3 
can be effectively reached. In fact, for atoms (and 
many other not so small systems), the total amount of 
matter is very small and leads to t^ v rj^. However, 
if /' is sensitive to ultra-low densities such as p M , the 
physical metric <? p „ = r/^/f departs from ij^ when 
the matter density m e \ip(x)\ 2 drops below p^. This 
is at the heart of the matter instabilities found in llUlPJ. 



On the other hand, the dependence of g M „ on the local 
matter-energy distribution also imposes constraints on 
the sources. Since the Ricmann tensor involves second- 
order derivatives of the metric, it follows that a smooth 
geometry must be one with continuous and differentiable 
/' up to second order, d 2 g iiv ~ (<9 2 /')V ~ (df') 2 t^, 
which implies that T must be differentiable up to that 
order at least [17]. An immediate consequence of this 
restriction is that the naive spherically symmetric con- 
structions of GR with matter profile of the form 



Po if r < R 
if r > R 



(8) 



are, strictly speaking, not valid for a generic Palatini 
f(R). This idealized construction, used as such in [H, 
can be easily adapted to Palatini f(R) by conveniently 
modifying the density profile near the surface to smoothly 
interpolate between the interior po and the exterior vac- 
uum. Otherwise, the discontinuous density profile would 
yield undesired singularities. Nonetheless, intuition tells 
us that such a refinement should be (from a practical 
point of view) unnecessary for theories with high curva- 
ture corrections, for which the contribution to the metric 
due to the matter contained in the outermost layers of the 
star should be negligible and then one could effectively 



use GR. However, it is precisely the existence of singu- 
larities on the surface of stellar objects that motivates 
this paper. In the next sections we will study in detail 
the properties and nature of such (highly unexpected and 
counterintuitive) singularities . 



III. STELLAR STRUCTURE EQUATIONS 



Using the following line element 



ds 2 = -A(r)e 2 ^dt 2 + 



1 



A(r 



■dr 2 



(9) 



with A(r) = 1 — 2M(r)/r, the combinations G\ ± G r r of 



lead to 

2 
r 



r 



Mr 



K 2 (P + P) 

f'A(r) 
k 2 ( p + V/2) 



r 
r 



+ A(r) 



f 

•I rr 

T 



(11) 



f'r [r-M(r)] 
/' r[r - 2M(r)] 



where /' = df/dR, f' r = df'/dr, f'„ = d 2 f/dr 2 , W = 
dip/dr, M r = dM/dr, and V = TZf - f. These two 
equations, plus the conservation equation 



dP , , ( , A r 



(12) 



plus an equation of state p = p(P) suffice to determine 
the structure of static, spherically symmetric objects. 
The equations corresponding to GR are recovered in the 
limit /' = 1. 

In order to proceed with the algebraic manipulations 
and extract information about the new physics hidden in 
equations (|10[) . (fTTj) , and p2[) . we must substitute (jTT 
and (JTTJ) into (T5]) to get 



dP 
dr 



(p + P) 



(n 2 P - V/2) M(r) 



A(r) [ /'A 
A{r)f r (lK i r-M(r) 
,4/' 



- (13) 



A /' 



2M(r)} 



where A = (2/r + f' r / '/')• The terms in the second line 
of this equation contain the main contribution due to the 
modified lagrangian. To deal with them, it is useful to 
express the matter terms f' r and f rr in terms of P r and 
P rr explicitly. Since f r = f"(dK/dT){dT/dP)P r , we 
can use ^ to find 



m 

dT 



Kf" - r 

This way we can express f' r as f' r = fpP r , where 

n 2 f" f dp 
Kf" -f\ dP 



(14) 



(15) 
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and we have used T = 3P — p. Similarly, one finds that 
f' rr can be expressed as f' rr = f' P P rr + f PP P 2 , where 



f'f 



K 4 fr 



PP 



dp 



w-rfv dpj (TZf-f')dP 2 

(16) 

Fortunately, f' rr does not appear in (JT3J) and the depen- 
dence on f r leads to a quadratic algebraic equation on 
P r which can be solved inmediately. The result can be 
expressed as follows: 



dP 

dr 



P 



(0) 



{I -a) 



where we have defined 
p(0) ^ (P + P) 



r[r-2M(r)] 
_ (P + P)f'p 



M{r) 



1+PP, 



K 2 P 



(0) 



(17) 



2j 2T 



f 



3r 2 [r-2M(r)} ( f' t 



2(1 -a) 2 



(18) 
(19) 
(20) 



Note that Pr°^ coincides, up to the factor V/2, with the 
result corresponding to GR. The terms a and [3 are pro- 
portional to f'p/f and {f'p/ f') 2 , respectively, and vanish 
in the limit of GR, /' = 1. Using (|17[) it is now possible 
to compute P rri which will be needed in (fTTj|) and (fTTjl . 
In terms of P r and P rr , those equations can be cast as 



1 



2/' 7 



n 2 {P + p) 



2 r 



P r M r 




We now have all the equations needed to study the struc- 
ture of spherically symmetric, static configurations in 
Palatini f{R) theories written in explicit form. 



IV. POLYTROPES IN f(R) = R±\R 2 

We will now study the properties of equations (TTT)) . 
([2"Tj) . and (|22[) for a matter distribution with polytropic 
equation of state 



.<'■' = I* 



l/ 7 



7-1 



(23) 



in the region close to the surface of such an object, i.e., in 
the limit P — > 0. Our gravity lagrangian is characterized 
by the Planck length scale A = Ip, which is related to 
a density scale p\ = (k 2 A) _1 ~ 2 • 10 92 g/cm 3 (we have 
used k 2 — 8nG/c 2 ). For these models, we find 



1 ± 2Ai? = 1 T 2- 



f'p = T 



2(3 -p P ) 



f'f 



pp 



where T = 3P — p and 
1 



Px 
2ppp 

Px 



1 



Pp 



Ppp 



P 



7 — 1 \ K 



(1-7) (P 



(l-7)/7 



K 



(l-2 7 )/ 7 



(24) 
(25) 
(26) 

(27) 
(28) 



Since we are interested in the regions near the surface, 
we now study the leading order of the different quanti- 
ties appearing in (JTTJ), (f2"Tj) , and (j2"2"|) . For very low pres- 



sures, we find that p m 



p \!/7 



^ 1 /P\(!-7)/7 

Pp ~ W \k) 
and ppp = fc^ll (jf) 7 ^ 7 remains unchanged. Us- 



ing this, p4]l. (|25j). and (|26J) become 



f'p 



1±2 



P 



P\ 
2 1 



1 ± 



P 



P 

Px \K 

(l-7)/7 



1/7 



f'f 



px \K 



2 (l- 7 ) (P 



pp 



(l-2 7 )/ 7 



px (jK) 2 \K / 
We also find that fl8|). (fT9|) . and ([20]) become 

1/7 



p(0) 



r{r - 2M tot ) 



1 1 



P 



K 

(2-7)/7 



pxjK \K 
Qr 2 [r-2M{r)\ 



1 



(1 - a) 2 {jK Px f 



(3-2 7 )/ 7 



(29) 
(30) 
(31) 

(32) 
(33) 
(34) 



With these results it is easy to check that, when P — > 
and for 7 < 2, a -> and /3P r (0) -> 0, and (pU) boils 

down to P r ps —Pr°^ oc (^) 1 ^ 7 , which is well-behaved 
and vanishes in this limit. Knowing the behavior of P r , 
it is easy to verify that f P P r , P r r, and f'pP rr are well- 
behaved for 7 < 2. However, the term f PP P 2 of (f2Tj) 
and (|2"2")) is problematic [lj|. In fact, from IpTF]) and 
it follows that 



f'ppPr 



2 (I-7) 



M? ot 



p x ^K) 2 r 2 {r-2M tot ) 2 



(3-2 7 )/ 7 



(35) 
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This term diverges for 7 > 3/2 when P — > 0. The 
existence of this divergence was first reported in Q, 
where the matching conditions between the internal 
and the external (Schwarzschild) metric were analyzed. 
Since this term is contained in ip r and M r , and such 
terms appear in the definition of the Riemann tensor, it 
was concluded that the geometry becomes singular near 
the surface of polytropes with index 2 < 7 < 3/2. Since 
the physically interesting case 7 = 5/3 (degenerate, 
non-relativistic fermion gas) lies within this interval, 
this result was regarded as a serious theoretical concern 
about the viability of Palatini f(R) theories. We will 
study now whether the (mathematical) divergence found 
in (f35|) is physical or not and how stringent is its effect 
on the viability of Palatini f(R) theories. 

First of all, for the models considered here, the diver- 
gence of f' PP P 2 is entirely due to the contribution coming 
from the polytropic equation of state. This can be seen 
from the definition of f' pp in (JIBJ) [or from (J^BJ) or ([51])] 
since dp/dP and d 2 p/dP 2 diverge as P ~ + irrespec- 
tive of the function f(R). If dp/dP and d 2 p/dP 2 were 
finite as P — > 0, the divergence would not arise. Second, 
a glance at the dimensionless (and exact) contribution 
^-{p+P) 2 of f' PP P 2 indicates that this term is strongly 

suppressed everywhere by the factor l/p\ ~ 10~ 92 . (This 
also happens with all other correcting terms, which in- 
dicates that for our Planck scale corrected theory the 
interior structure of stars is virtually the same as in GR. 
In [Tol | , however, the value of A was chosen in such a way 
that the interior of compact stars was indeed affected by 
the full f(R) dynamics.) The exception occurs strictly in 
the limit P — > 0. Re-expressing that dimensionless factor 
in terms of the density, we find 

^(P + P)^^$P^. 06) 

Since for 7 > 3/2 this term diverges as p — * 0, it is worth 
looking at the value at which it becomes of order unity. 
This happens at p — p s , where 

Ps=[—) ■ (37) 

For non-relativistic neutrons, p <C 6 • 10 15 g/cm 3 , 
the ideal Fermi gas approximation leads to a poly- 
tropic equation of state of index 7 = 5/3 and K — 

(3 2 / 3 ir 4 / 3 /5){h 2 /m S n /3 ) « 5 • 10 9 (in c.g.s. units), where 
m n is the neutron mass. The resulting p s is on the or- 
der of p s ~ 10~ 210 g/cm 3 [l9(. This density is, by far, 
well below any physical density one can imagine. In fact, 
for a free electron whose wave function is spread over 
the entire universe, the ratio m e /R 3 Jniv is on the order 
of ~ 10~ 118 g/cm 3 . This means that a single electron 
outside of this idealized polytrope is more than enough 
to cure this singularity (of course, we assume the elec- 
tron wavefunction uniformly spread over a spherical shell 



to respect the symmetry of the problem) since its mere 
presence rises the average matter density 92 orders of 
magnitude above the critical scale p s . Therefore, the ex- 
istence of a curvature singularity at such extremely low 
densities should be regarded as unphysical, as an arti- 
fact of the idealized equation of state used. One should 
have in mind that an accurate EoS at laboratory densities 
(let alone at the density scale found here) is very com- 
plicated to derive, because electrostatic interactions and 
other subtle effects mask the simpler statistical effects of 
the idealized Fermi gas approximation [l3|. The poly- 
tropic EoS should therefore be used within its expected 
regime of validity. 

V. DISCUSSION AND CONCLUSIONS 

This work was motivated by an extremely interesting 
and puzzling result: the existence of curvature singulari- 
ties in regions of very low curvature affecting all families 
of f(R) lagrangians in Palatini formalism. If such sin- 
gularities were physical, all f(R) lagrangians except that 
of Hilbert- Einstein should be ruled out. And that would 
be so for arbitrarily small departures from f(R) = R. 
Such a result would have had very strong physical conse- 
quences, since it would have singled out GR as the priv- 
ileged f{R) lagrangian. Our analysis, however, shows 
that the (mathematical) low density/curvature singular- 
ities found in Q depend intimately on the particular EoS 
chosen and have different strength depending on the f(R) 
model chosen. In the model that we studied, character- 
ized by the Planck length squared, a single electron would 
suffice to cure all stellar singularities in the Universe. 
If instead of the Planck length we had chosen a larger 
length scale (smaller density scale p\), the surface sin- 
gularity could occur at higher densities, perhaps within 
the regime of validity of the polytropic EoS, and could 
not be removed by such simple means. In this sense, the 
doubts raised in [3] and [l(| about the viability of Pala- 
tini f(R) models are well justified. In those works, the 
value of A was chosen on phenomenological grounds to 
be A = (0.15 km) 2 , which leads to p\ ~ 2 • 10 18 g/cm 3 . 
This density scale is enough to pass all solar system weak 
field tests [14| since /' « 1 to high accuracy everywhere. 
However, it leads to p s as 4 • 10 12 g/cm 3 , which indicates 
that the divergent terms of M r and tp r begin to grow well 
within the region of validity of the polytropic EoS. The 
study of polytropes carried out here and in Q and [Io| 
can thus be used to place bounds on the parameter A 
of the model under consideration and also in other mod- 
els. If, for instance, we assume that the polytropic EoS 
should not be trusted below p s ~ 10~™ g/cm 3 , we find 
that 

A < IO 4 -' 1 / 3 cm 2 , (38) 

which places explicit constraints on the allowed values 
for A. 
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To conclude, the important point was to show that 
not all lagrangians are ruled out by this theoretical 
experiment with polytropes. We have confirmed that the 
surface singularities exist, but their physical character 
depends on the details of the model and, as expected, 
ultraviolet corrections are allowed for suitable choices 
of parameters. In particular, Planck scale corrected 
models are save from such singularities. This is of 
great importance to the light of the recently discovered 
relation existing between non-perturbative approaches 
to quantum gravity and Palatini f(R) theories [HI]. It 
turns out that the only consistent homogeneous and 
isotropic cosmology that can be constructed using the 
techniques of Loop Quantum Gravity and which is 
free of the Big Bang singularity [l(| can be put into 
correspondence with a Palatini f(R) lagrangian with 



ultraviolet, Planck scale, corrections. This provides a 
fundamental and so far non-existing theoretical justifi- 
cation for Palatini f(R) theories, which have only been 
studied on phenomenological grounds in relation with 
the cosmic speedup. We hope that this new avenue 
of research helps us better understand the properties 
of quantum gravity and space-time near cosmological 
and black hole singularities using the more familiar 
techniques of modified gravities. 

The author thanks E.Barausse, T.P. Sotiriou, and 
J.C. Miller for their critical reading of this manuscript 
and clarifications on their works. Special thanks go to 
F.Barbero for very useful comments and discussions on 
this subject. This work has been supported by MICINN 
through a Juan de la Cierva postdoctoral contract. 
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